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Abstract. In this article, I review the most recent results in radiative penguin decays from the B 
factories Belle and BaBar. Most notably, I will talk about the recent new observations in the decays 
B — » (p/uj)j, a new analysis technique in b — > sy, and first measurements of radiative penguin 
decays in the B® meson system. Finally, I will summarize the current status and future prospects 
of radiative penguin B physics at the B factories. 
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1 Introduction 

Now, at the verge of the Large Hadron Collider (LHC) 
becoming operational, an exciting new era of discovery 
in high-energy particle physics is anticipated. There is 
a whole series of new physics models attempting to 
fix some of the shortcomings of the Standard Model 
(SM) of elementary particle physics, most notably the 
hierarchy problem of the Higgs sector. Low energy su- 
persymmetry (SUSY) is probably the most popular 
extension of the SM. But all these new models intro- 
duce new free parameters. Current precision experi- 
ments can be used to constrain the allowed space of 
these additional new parameters and therefor guide 
the searches at the LHC. 

One very interesting class of decays are so-called 
penguin decays where the leading order SM decay con- 
tribution is described by a one-loop Feynman diagram. 
Here, I will concentrate on penguin decays of B mesons 
where a photon is radiated off. 

Unknown particles that do not exist in the SM can 
contribute considerably to these decays, e.g. by ap- 
pearing in this loop, and thus shift measurable observ- 
ables. But if the measurements agree with SM based 
theoretical predictions, constraints on various beyond 
the SM physics scenarios can be extracted. 

There is a large interest in B meson decays since 
they provide an excellent window into the physics of 
quark mixing. In the SM, this is parameterized by four 
numbers that appear as the four independent parame- 
ters of the Cabibbo-Kobayashi-Maskawa (CKM) quark 
mixing matrix. A stringent test of the SM is performed 
by overconstraining these four parameters with inde- 
pendent measurements, of which many can be done by 
studying B meson decays. 

Hundreds of millions of B mesons have been pro- 
duced in the process e + e~ — > T(4S) — > BB at the two 

a Email: kkoeneke@alum.mit.edu 



energy-asymmetric e + e~ colliders KEKB and PEP- 
II. The measurements reviewed in this paper are per- 
formed with the Belle and BaBar detectors p][2] situ- 
ated at these two B factories. In Sec. [2j I review com- 
mon reconstruction techniques. In Sec. [3l I will talk 
about the decays B — > {p/co) 7. I introduce a new ex- 
perimental approach to measuring the decay b — ► s^Q 
in Sec.[4j In Sec. [5) I report on results of radiative pen- 
guin Bg decays using a data sample collected at the 
T(5S) resonance. Finally, I conclude in Sec. [HI 



2 Common reconstruction techniques 

B mesonfl are reconstructed in two standard vari- 
ables, rriEs & n d AE. The first is defined as titles — 

\J Z ~ IPb| 2 ; where s is the center-of-momentum (CM) 
energy of the e + e~ collision and p* B is the measured 
three- vector of the reconstructed B meson, also in the 
CM frame. The second important variable is defined 
as AE = E* B — y/s/2, where E* B is the measured CM 
energy of the reconstructed B meson. 

Due to the small mass difference between m r ( 4S ) — 
2itib, the two B mesons are almost at rest in the CM 
frame. They decay thus rather spherically symmetric 
in that frame. This is different for the so-called light- 
quark continuum background (e + e~ — * //, where / = 
u, d, s, c, t) where the light quarks (or leptons) have a 
large momentum in the CM frame. Their hadroniza- 
tion and decay is thus relatively collimated along the 
decay axis. These event shape differences are used in 
various ways in all analysis described below. 



Charge conjugation is implied throughout this article. 
2 All particles mentioned in the experimental part of this 
paper are actually "candidates" in the sense that their true 
identity is only assumed. 
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3 The Unitarity Triangle and B — ► (p/oj)^ 

Besides the generic interests in radiative penguin de- 
cays mentioned in the introduction, this decay is also 
interesting due to its implications to measuring CKM 
matrix elements. If the experimental branching frac- 
tion of B — > {p/lo) 7 is known, the ratio with that of 
B — > K*j is proportional to |Vtd/^ts| 



g(B->p 7 ) _ q 

B(B^K'~f) °P 



Vfd 



(l-mg/M 2 ) 3 
C 1 - 



,/M 2 



(1) 



[1 + AR{p/K* 



Table 1. The branching fraction (B) for each B — > (p/uj) 7 
mode is shown for the results from BABAR [10] and the new 
preliminary results from Belle QT], including their signifi- 
cance (£) in standard deviations (systematics included). 



Mode 


B ± cr stat ± t 
BABAR 


r sys (10- 6 ) (S) 
Belle 


p+7 
P°7 


I.IO+0 33 ± - 09 ( 3 - 8 ) 

0.79lo;2o ± 0.06 (4.9) 
0.401^20 ± °- 05 ( 2 - 2 ) 


0-°6_o. 28-0.08 \ 6 - Z ) 

0.76 ±0.17 ±0.06 (4.9) 
0.42^;^ ± 0.04 (2.6) 


(/»/w)7 
PI 


1.251^24 ± °- 09 ( 6 - 4 ) 
1.36+0 2? ± 0.10 (6.0) 


1.13 ±0.20 ±0.11 (5.9) 
1.19 ±0.24 ±0.12 (5.5) 



Here, £ = £±(0)/£± (0) is the ratio of form factors 
computed in Heavy Quark Effective Theory (HQET), 
S p = 1(1/2) are isospin weights for the charged (neu- 
tral) p meson and AR(p/ K*) is a dynamical function 
calculated for example in |4J which accounts for dif- 
ferent dynamics in the three decays like vertex, hard- 
spectator, weak annihilation or exchange contributions. 

A recent calculation based on Light-Cone Sum Rules 
determines £ p = 1.17 ± 0.09 and = 1.30 ± 0.10 [5]. 
AR(p/K*) is usually computed to be in the neigh- 
borhood of 0.1 ± 0.1. A new calculation that also ac- 
counts for the additional W annihilation diagram of 
the B + — » p + j decay is also available [BJ. 

The same ratio of CKM matrix elements is also 
measured by a physically very different process, B mix- 
ing. The frequency of the oscillation of a B meson turn- 
ing into a B meson is related to the mass difference of 
the light and heavy B meson eigenstates Am. The ra- 
tio of these mass differences between the B% and the 
B° s is proportional to |Vt<j/Vt s | 2 . The B factories have 
measured Arrid very precisely [7] and Am s has recently 
been measured by the CDF collaboration [5]. 

The first observation of B — > (p/oj) 7 decays was 
reported by Belle regarding the B° — ► p°7 channel 
and the combined B — > (p/cu)j channel [5]. BaBar 
also reported recently the observation of the combined 
B — > (p/ui) 7 and B — > pq modes and also evidence for 
B° — > p°7 and first evidence for B + — > p+7. This re- 
sult is based on a dataset of 347 million BB events [10] . 

The high-energy photon is enforced to not originate 
from a tt° or r\ decay by means of a two-dimensional 
likelihood comprised of the energy of any other pho- 
ton found in the event and the two-photon invariant 
mass. Continuum background is rejected with a high- 
dimensional neural network, optimized for each mode 
individually. These neural networks combine event shape 
information, separation of the decay vertices of the two 
B mesons in the event along the z axis, decay angular 
information, information related to tagging of the fla- 
vor of the other B meson in the event, and other quan- 
tities. The signal extraction is performed with an un- 
binned maximum likelihood fit in four (five) dimension 
for the individual B — > pj (B — > 077) channels, where 
the dimensions used are wies, AE, a transformation of 
the neural network output, and the cosine of the helic- 
ity angle of the vector meson (and the sine of the sec- 
ond independent angle in the three body u — > 7r+7r~7r° 



decay). For the extraction of the combined results a si- 
multaneous fit is used with the additional constraint 
on the decay widths / n g_ >p + 7 = 2Pg^ ( o°7 = 27^B_ >a , 7 . 

The resulting branching fractions are listed in Ta- 
ble [1] While writing these conference proceedings, the 
Belle collaboration presented an updated preliminary 
analysis of these decays at the Lepton-Photon 2007 
conference performed with a data sample of 657 mil- 
lion BB events [TT]. Evidence is found for the two pj 
modes and both combined modes are observed, but 
the B° — > ui-f mode still eludes detection. These pre- 
liminary branching fractions are also listed in Tabic [T] 

Combining both measurements to a new prelim- 
inary world average and using the measured results 
of the B — > if*7 decay [7], which includes the lat- 
est BaBar measurement [H] , and applying Eq. Q] with 
input from [5][6] yields the results 

\Vu/V ts Zi = 0.194lHl!(exp.) ± 0.014(th.) 
\Vtd/V t s\p = 0.201 ± 0.016(exp.) ± 0.015(th.) 

for the p/tu mode and p+ / p° mode, respectively. 

Both of these results are in excellent agreement 
with the result extracted from B^JB® oscillations re- 
cently reported by the CDF collaboration [5] 

\Vt d /V tg \TZ/A ms - 0.208t™°i(expOtS(th.)- 

The implications of the |Vtd/Vt s | measurement using 
the two p modes only is shown in Fig. [TJ compared 
with the CDF result. It is a remarkable success of the 
Standard Model that these two very different physics 
processes agree that beautifully. 



4 A new approach to b — > 37 

Besides the possibility of constraining extensions of the 
SM, this decay is interesting because the high-energy 
photon carries information about the b quark inside 
the B meson. If we consider for a moment the decay 
of a free b quark to a free s quark and the high-energy 
photon, then this photon has a specific energy since 
this would be a two-body decay. As quarks are not 
free, the energy of the photon is smeared. But still, 
its first moment is related to the mass of the b quark 
and its second moment to the Fermi motion of the b 
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Fig. 1. Constraints on the far side of the unitarity trian- 
gle extracted from the B — > py measurements reported in 
this paper (green) compared to B oscillation measurements 
(yellow), computed with the UTFit code [T3lfT4] . 

quark inside the B meson. Extracting these quantities 
from the high-energy photon energy spectrum feeds 
into the determination of the CKM matrix element V u b 
since the extracted parameters are universal for the B 
meson and are needed for the determination of this 
CKM matrix element from semi-leptonic B decays. 

The three recently published theoretical SM branch- 
ing fraction predictions at NNLO [T5llT6}|T7] show the 
big interest of the physics community in this decay. 
These precise calculations (up to 7% uncertainty) al- 
low for a very good comparison with the measurement 
and thus a stringent test of the SM can be performed. 

Traditionally there are two experimental approaches 
to this measurement, a fully inclusive and a sum-of- 
exclusive. In the fully inclusive approach only the high- 
energy photon of the signal B decay is reconstructed 
and the event is tagged as a T(AS) — > BB event by 
requiring the existence of a lepton (e or /j) coming 
from a semi-leptonic decay of the other B meson (the 
"tag B") in the event. This method is theoretically 
very clean but has two experimental drawbacks: high 
background and the measurement of the high-energy 
photon energy spectrum in the T(4S) rest frame in- 
stead of the B rest frame. The sum-of-exclusive ap- 
proach reconstructs the decay in as many final state 
modes of the s quark fragmentation as possible. Thus 
the photon energy is measured in the B rest frame and 
the background is much reduced w.r.t. the fully inclu- 
sive method. However the theoretical interpretation is 
more difficult due to the uncertainty of the missing 
fraction of the s quark hadronization. 

A new experimental approach presented by BaBar 
is combining the advantages of both analysis. The tag 
B meson is fully reconstructed in an all-hadronic final 
state, e.g. B — > Dir. Only the photon is reconstructed 
from the signal side B meson. Due to the knowledge 
of the initial beam conditions and the full rcconstruc- 



Table 2. Measurements of the first and second moment 
of the photon energy spectrum in the 6 — » sy decay, de- 
pending on the low cut on the photon energy _E^ ut . 



J-iCUt 

t 1 


First moment 


Second moment 


(GeV) 


Value ±<7 sta t ± CTsys 


Value ±(7 s tat ± OSys 


1.9 


2.289 ±0.058 ±0.026 


0.0334 ±0.0124 ±0.0065 


2.0 


2.315 ±0.036 ±0.020 


0.0265 ±0.0057 ±0.0024 


2.1 


2.371 ±0.025 ±0.011 


0.0142 ±0.0037 ±0.0013 


2.2 


2.398 ±0.016 ±0.006 


0.0092 ±0.0015 ±0.0010 


2.3 


2.427 ±0.010 ±0.007 


0.0059 ±0.0007 ±0.0003 



tion of the tag B meson, the four momentum of the 
signal side B meson is known and thus, the photon en- 
ergy can be measured in the B meson rest frame. The 
fully tagged event also reduces the background consid- 
erably. No uncertainty from the s quark hadronization 
is present since only the high-energy photon is recon- 
structed and not the hadronization products of the sig- 
nal side s quark. The drawback of this new method is 
the small reconstruction efficiency of only about 0.3%. 

After the tag B meson reconstruction, the signal 
side high-energy photon is required to not be a prod- 
uct of a 7T°, 77, or p decay. The continuum background 
is suppressed by means of a Fisher discriminant which 
is built with 12 inputs, mostly based on event shapes. 
The remaining continuum background does not peak 
in tiies and can thus_ be subtracted. The left over 
e + e~ — > T(4S) — > BB background is mostly due to 
a 7T° or 77 decay where the signal high-energy photon 
is not identified as a decay product of these particles. 
The signal yield is determined from an individual fit 
to ities for each 100 MeV photon energy bin. 

The branching fraction for a minimum energy of 
the photon of 1.9 GeV is determined on a BaBar dataset 
with an integrated luminosity of 210 fb _1 as 

B{b^s~f) [E n > 1.9 GeV] = (3.66 ± 0.85 ± 0.59) xl0~ 4 . 

Extrapolating this to the theoretically preferred photon- 
energy cut of 1.6 GeV using [I8j gives 

B(b^> s 7 ) [E y > 1.6 GeV] = (3.91 ± 0.91 ± 0.63) xl0~ 4 . 

The systematic uncertainties can be reduced with higher 
statistics. The resulting first and and second moments 
for different photon energy cuts are shown in Table [21 
With this new preliminary BaBar measurement, a 
third independent approach to measure b — > sy is ex- 
plored. Since the data samples of the three measure- 
ments are independent, they can be combined to one 
result. Currently, this new method does not reach the 
precision of the other two. But I expect that unlike the 
other two, this approach will benefit from more statis- 
tics, especially from very large datasets in the multi 
inverse attobarn era expected at super-f? factories. 

5 Results from the T(5S) run 

KEKB ran also at a higher CM energy equivalent to 
the mass of the T(5S) meson. The Belle detector col- 
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lected 23.6 fb at this configuration, corresponding 
to about 2.6 million B® mesons produced. The number 
of produced B® mesons per inverse femtobarn is small 
for two reasons. First the production cross-section of 
the T(5S) resonance is with 0.3 nb about three times 
smaller as that of the T(4S) resonance. And only about 

20% of the T(5S) mesons decay into B® pairs. The 
majority of strange B mesons are actually created as 
B° s * mesons and they decay via B®* — > B® 7 with a soft 
photon. Besides this soft photon, B® reconstruction is 
performed as in the T(AS) — ► BB case. The photon is 
too soft to be reliably reconstructed, so it is ignored 
and only visible in a ss 50MeV downward shift in AE. 

5.1 B° s -> 07 

This decay mode is the brother of B° — * if* 7 , just the 
spectator quark is different. The SM theory prediction 
for the branching fraction is (4 ± 1) x 10~ 5 [6]. A recent 
NNLO calculation gives (4.3 ± 1.4) x 10" 5 [19]. 

Two oppositely charged kaons are reconstructed 
and combined to form the <f> candidate. A cut on the 
invariant mass of the two kaons is applied with a 2.5a 
window around the nominal <p mass. The light quark 
continuum background is suppressed utilizing event 
shape information in the form of modified Fox- Wolfram 
moments. The photon is required to not originate from 
a 7T° or an r\ decay. The final signal yield is extracted 
with a three-dimensional unbinned maximum likeli- 
hood fit. The three variables entering in this fit are 
rriESi AE, and the cosine of the helicity angle of the 
<j) meson. 18 signal events are found which yields 

B(5 s °^0 7 ) = (5.7l^; 2 7 )xlO- 5 . (6) 

The significance of this result, including systematic un- 
certainties, is 5.5cr. This is the first observation of a 
radiative penguin decay of the B® meson. 

5.2 B° s -> 77 

This decay proceeds in the SM at leading order via 
a penguin fusion diagram between the b quark and 
the s quark. It is like the Feynman diagram of the 
b — > S7 decay turned by 90° with an additional photon 
attached. The SM theory predictions are in the range 
(0.5 - 1.0) x 10" 6 [201I2T] . This low branching fraction 
can be enhanced by up to one order of magnitude in 
several models extending the SM, e.g. in 4 th quark 
generation models |22j or in supersymmetric models 
with broken R-parity [23]. No discovery of this decay 
has been reported as of today. 

The reconstruction is performed similar to the above 
case. Especially the two photons are required to not 
be consistent with the decay of a ir° or rj meson. The 
signal is extracted with a two-dimensional unbinned 
maximum likelihood fit where the two variables are 
rriES and AE. No signal is seen and an upper limit is 
set at a 90% confidence level of 

B (B° s -► 77) < 8.6 x 10~ 6 at 90%. (7) 

This is as of today the best limit for this decay mode. 



6 Conclusions 

New measurements of B — > (p/ui)j are available. The 
agreement with B® oscillation measurements is excel- 
lent and a non-trivial test of the SM succeeds. The ex- 
perimental precision will improve with larger datasets. 
A new reconstruction method for the b — ► sj decay is 
presented. This method will fully benefit from datasets 
from super B factories unlike previous reconstruction 
methods. The first radiative penguin B decay of the 
B® meson, B® — > <jrj, has been observed by Belle using 
data collected at the T(5S) resonance. All measure- 
ments presented have a large impact on SM tests and 
constraints of beyond the SM physics. The rich field 
of radiative penguin B decays will provide more input 
to constraining the physics beyond the SM. 
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